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, Stochastic Schrodinger equations for quantum trajectories offer an alternative and sometimes 

superior approach to the study of open quantum system dynamics. Here we show that recently 
established convolutionless non-Markovian stochastic Schrodinger equations may serve as a pow- 
CN| , erful tool for the derivation of convolutionless master equations for non-Markovian open quantum 

systems. The most interesting example is quantum Brownian motion (QBM) of a harmonic os- 
^ ' cillator coupled to a heat bath of oscillators, one of the most-employed exactly soluble models of 

open system dynamics. We show explicitly how to establish the direct connection between the ex- 
act convolutionless master equation of QBM and the corresponding convolutionless exact stochastic 
Schrodinger equation. 
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I. INTRODUCTION 



Current attempts to push quantum coherent dynamics further towards the macroscopic Q, Q or towards systems 
residing on more and more particles [Sj have to overcome decohering influences of the environment. Typically, the 
system in question becomes entangled with environmental degrees of freedom and thus coherence is lost. Even if 
seemingly weak, couplings to an environment may have dramatic consequences for coherences. Thus, any theoretical 
Oh, modeling of such quantum dynamics should take into account environmental effects. 

+-> ' Traditionally, open quantum systems are described in terms of their reduced density operator and the dynamics 
thereof. Ideally, maps, propagators or even evolution equations are derived that account for environmental effects 

mum- 

Q-i The last decade has seen tremendous progres s in the development of stochastic Schrodinger equations to describe 
open quantum system dynamics jj, 0, 0, HU LUl 03 • I n this framework, the reduced density operator is obtained 
from an ensemble mean over the pure state solutions - "quantum trajectories" - of a stochastic Schrodinger equation. 
While known to exist in the case of Markovian open system dynamics (when the ensemble dynamics is governed by a 
Lindblad master equation |l4. Il5|). it is only fairly rec ently that quan tum trajectories were extended to cover more 
general non-Markovian open systems |H lmlllL lmlMl2l L l2^IH l2ll25t . Interestingly, in our approach [i^ ITl Il9| 
the non-Markovian stochastic Schrodinger equation is derived directly from an underlying microscopic system-plus- 
environment model without referring to the corresponding evolution of the reduced density operator. This turns out 
to be an efficient way to tackle non-Markovian quantum dynamics, shedding new light on the difficulties encountered 
when dealing with memory effects. 

Non-Markovian dynamics usually means past contributions to the current time evolution: memory effects typically 
enter through integrals over the past |26l l27|. Nevertheless, under certain circumstances, non-Markovian dynamics 
may be cast into a time-local, convolutionless form where the dynamics of the open system state is determined by 
the state at the current time t only [28l 12^. l3fj| . Then, non-Markovian effects are taken into account by certain time 
dependent coefficients that may replace the memory integrals usually encountered. 

Examples of this class are models of open quantum system dynamics that may be treated without any ap- 
proximation. The famous d amp ed harmonic oscillator bilinearly coupled to a bath of harmonic oscillators 
|28l I31L I32I l33l I34L l3ol l36l I3H l38| allows for a convolutionless treatment without any approximation. Such model 
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systems allow us to test approximation schemes and, very often, they enable us to pin down important qualitative 
features of the dynamics of more realistic open quantum systems clearly. 

The purpose of this paper is to show that our non-Markovian quantum trajectory approach, besides its powerful 
use in numerical simulations, provides an elegant way to derive non-Markovian convolutionlcss master equations for 
open quantum systems. We exemplify this idea by deriving the exact convolutionless non-Markovian master equation 
for the famous quantum Brownian motion model of a harmonic oscillat or 13 4 ). l37| . This QBM master equation has 
been extensively studied in different contexts [U El IH Eftllllll Ellilli^. Our approach is 
different from all the previous work. We show in this paper how to establish the direct link between our convolutionless 
stochastic Schrodinger equation and the known convolutionless master equation. The key ingredient for the derivation 
is a Heisenberg operator approach to our non-Markovian stochastic Schrodinger equation [24J . 

Our paper is organized as follows: in Section II we review the " non-Markovian quantum state diffusion " stochastic 
Schrodinger equation, sketch its derivation and display different approaches to a convolutionless formulation. Next 
we show how to derive master equations from the stochastic Schrodinger equation and give some simple examples. 
The main part of the paper follows in Sections IV and V where we discuss both the stochastic approach to the soluble 
QBM-model and the derivation of the convolutionless master equation based on the stochastic approach. Unavoidable 
technical details are left for two appendices. 



II. STOCHASTIC SCHRODINGER EQUATION 



Here we review the Markov and, more importantly, our non-Markovian stochastic approach to open quantum system 
dynamics |l8l [19| . The dynamics of the reduced density operator p t is obtained from a stochastic evolution equation 
for pure states tpt, such that upon taking the ensemble mean A4{. . .} over the pure states, the reduced density operator 
is recovered: 

Pt = M{\ilH)(ilh\}. (1) 



A. Usual Markov limit: Lindblad theory 

In the standard Markov case, often encountered in quantum optical a ppli cations, the evolution equation for the 
reduced density operator of the open system takes the Lindblad form [T^.ll5j 

dtPt = ~[H,pt] + ~ ([L, Pt tf] + [Lpt,tf]) . (2) 

We here restrict ourselves to the simplest case of a single non-unitary contribution involving a single operator L (the 
Lindbladian) . In general, a sum or integral over such Lindbladians may appear. Here and throughout the paper, we 
denote with H the Hamiltonian that generates the otherwise isolated, unitary evolution of the system. 

A stochastic Schrodinger equation "unraveling" Lindblad evolution @ is provided by the linear quantum state 
diffusion [Toll2^ | equation 

Otin =(-j; H + Lz t - \^L\ i> t (3) 

for unnormalized stochastic states ip t - The stochastic process z\ is complex white noise with zero mean (A4{z^} = 0) 
and correlations 

M{z t z*} = 5{t-s) and M {z t z s } = 0. (4) 

Equation is here understood as a stochastic equation in its Stratonovich sense, though differences to the Ito case 
only appear in their respective nonlinear versions, which are not relevant for this paper, see [lOj. As one may easily 
confirm, provided the quantum trajectories follow equation ©, the ensemble mean Q indeed evolves according to 
Lindblad's equation J2J. 



B. Non-Markovian stochastic Schrodinger equation 



In order to extend the Markov theory we use a standard model for open quantum systems in the system-plus- 
reservoir framework: a system with Hamiltonian H, coupled linearly to a large number of harmonic oscillators with 
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distributed eigenfrequencies u>\ and creation and annihilation operators b x ,bx- This total Hamiltonian can be written 



as 

H tot = H + H int + H hath = H + hJ2(9*xLb{ + g x ZSb x ) + £ hu x b{b x , (5) 

A A 

where L is the system operator providing the coupling to the environment, and g\ are coupling constants. 
The interaction between system and environment may be written in the more appealing form 

H iat = H(LB^ + L^B) (6) 

with a bath operator B consisting of contributions of all environmental oscillators, B — g\b\. It turns out to be 
convenient to change to an interaction representation with respect to the free bath evolution, such that instead of (JSJ 
we use 

Htot(t) = e* H ^ ht (H + tfinOe-^^ = H + % (LB\t) + l)B{t)) (7) 
for the total Hamiltonian with 

B(t) = e *^ b t Be -tH batll * = ^2g x bxe- iu ^. (8) 

A 

The dynamics of the open system is thus influenced by the bath operator B(t), whose statistical properties are 
captured in its correlation function 

/>OG 

a(t -s) = (B(t)B\s)) cnv = V | ff A| 2 e-^ ( *- s) = ft / duJ{w)e~^ t - s \ (9) 

a J ° 

here for a zero temperature bath (for finite temperature see below). The last equality in © defines the spectral 
density J(uj) of bath oscillators. 

If the dynamics of system and environment is such that the bath correlation function in @ may be replaced by a 
delta function, the reduced dynamics is Markovian and pt evolves according to the Lindblad master equation © . For 
a general correlation function a(t — s), however, memory effects of the environment may be important and matters 
become exceedingly more difficult. Such non- Markovian effects are known to be relevant in many situations, in 
particular at low temperatures and as soon as narrow energy splittings occur, as in tunneling processes [(|. They are 
also relevant for structured environments, as for instance the electro-magnetic vacuum in the presence of a photonic 
band gap material |44j |. Moreover, non-Markovian effects play a role in output coupling dynamics of atoms from a 
Bose Einstein condensate with the aim to build an atom laser . 

Also, it may happen that the dynamics appears to be Markovian, yet the evolution equation is not of the standard 
Lindblad class ©. This is the case, for instance, in the high-temperature limit of the standard quantum Brownian 
motion model to be discussed in Sect. IIVI In such cases, transient effects (initial slips, see also [46|,|42]]) are important 
and violate the semigroup property required for Lindblad-class evolution. 

The open quantum system described by the system-bath model JSj, allows to derive the linear non-Markovian 
generalization of ©, henceforth referred to as the "Non-Markovian quantum state diffusion" (NMQSD)- equation. 
The starting point is the formal solution l^t) of the Schrodinger equation for the total system 

*na t |* t ) = fl to t(*)|*t>, (io) 

in a special representation |18j . For simplicity, in this section we assume a zero temperature bath: the initial pure 
state of system and environment is 

|*o> = l^o) ® |0i) ® |0 a ) • • • ® |0a) ® • •• . (11) 

The system state |^>o) is arbitrary and all environmental oscillators start in their respective ground state 0a) - By 
using a Bargmann coherent state basis 0] for the environmental degrees of freedom: \z\) — exp{zAd A }|0A), and the 
resolution of the identity 

Ix= f^e~^ 2 \zx){zxl (12) 
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the total state |\&t) in (|10|) can be expanded as 

I**)- / ^e-l z l> t (z*))®|z), (13) 

where \z) — \z\) ® \z 2 ) ® ■ • • <S> \z\) ■ ■ •, d 2 z — d 2 zid 2 z 2 ■ ■ ■ d 2 z\ ■ ■ ■ and \z\ 2 — J2\ \ z \\ 2 - It follows from (|T3)l that the 
reduced density operator of the open system may be obtained by an ensemble mean over the system states \ipt(z*)), 

Pt = TVbathPt><M> t |] - J ^e~^ 2 \Mz*))(Mz*)\ = M{\iP t )(A\}- (14) 

The last equality defines the ensemble mean M{. . .} as the Gaussian integral J ^e~' z ' {. . .} over all coherent state 
labels z\. Using the Schrodinger equation (fTUf) . the pure states \ipt(z*)) of the system in (jT3|l were shown to satisfy 
the non-Markovian QSD equation (17L Il8| 

dtipt = ~H^t + Lz*^ - it / ds a(t - s)P^. (15) 

n j dz* 

where a(t — s) is the bath correlation function © and z* t = —i g^z\e lUxt is colored, complex Gaussian noise with 
M{z t } = M{z t z s } = and M{z t z*} = a(t - s). 

The linear non-Markovian QSD equation (|15f) was here derived for a zero-temperature heat bath. However, it can 
be shown that equations JQl and (|15|) remain valid also at finite temperatures provided we deal with Hermitian Lind- 
bladi ans, L = . In general, the finite temperature non-Markovian QSD equation takes a slightly more complicated 

form mm. 



C. Convolutionless stochastic Schrodinger equation 



Memory effects of non-Markovian evolution clearly make themselves felt through the integral over the past in 
equation (|15f) . involving the bath correlation function and a functional derivative of the current state ip t {z*) with 
respect to earlier noise z*. In many relevant cases, it is possible to replace that functional derivative by some time 
dependent operator O, 



Sz: 



0(t,s,z*)Mz*), 



(16) 



acting in the Hilbert space of the open system on the current state ipt(z*) |T^. We indicate that O may depend on the 
times t and s, and possibly on the (entire history of the) stochastic process z t *. Relevant examples of this replacement 
will be given shortly. One way to determine the operator 0(t, s, z*) in actual applications 0] is to insert Ansatz ll(j|) 
in 1|15|) and use the "consistency condition" 



' 6z* 



Sz: 



-dti>t{z*) 



together with the initial condition: 



for all s. 

One finds the formal evolution equation 



h 



d t O{t,s,z*) -- 
with the time-integrated operator 0(t, z*)\ 

0(t,z*) 



0(t = s,s,z*) = L 



H + Lz* t - tfO(t, z*), 0(t, s, z*) 



it 



SO(t, z*) 
5z* 



a(t — s)0{t, s, z*)ds. 



(17) 



(18) 



(19) 



(20) 



Equation Ijl9|l has to be solved with initial condition Ijl8(l for all s. As has been shown in previous publications [Isj . 
one can determine the O-operator in (|16(l for many interesting and relevant physical models. Moreover, approximate 
O-operators can always be found systematically and be easily implemented in numerical simulations |22l l'23l l53| . 
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Once the replacement (|16|) of the functional derivative by an operator is known - sometimes only approximately - 
the NMQSD equation (|15fl takes the more useful convolutionless form [Tj| 

dtM**) = + Lz* t - tfd(t, z*)^j i> t {z*), (21) 

where we used the notation 0(t, z*) from (121 )|) . see also The determination of ipt{z*) is now reduced to solving 
the simple time-local stochastic Schrodinger equation (|21|l . For a memory-less bath with a(t — s) = 7<5(i — s), we see 
from i|18|) and (|20|l that (5(i, 2*) = Sj-L holds. Then equation (|21|l reduces indeed to the Markov linear quantum state 
diffusion equation ©. 

We recall that (|21|) does not preserve the norm of the states ip t {z*). Thus, if only a numerical solution of (|21|l is 
possible, it is most often more advisable to use its nonlinear, norm-preserving version [Isj . 



D. Heisenberg approach 

There is an alternative approach to the convolutionless NMQSD stochastic Schrodinger equation l|2ip. see 0|. We 
find it most appropriate to work in terms of a stochastic propagator Gt{z*) for the states ipt{z*)- By definition, we 
have 

\M**)) = Gt(z*)\ik). (22) 

We use the unitary propagator U t for the total state (| 1 31) . satisfying 

ihd t U t = H tot {t)U t and U = I, (23) 

with H tot (t) from Q. Since \^ t ) = U t \^o) an d \ipt{z*)) = (z\^ t ) according to lJT5)l. the stochastic propagator G t (z*) 
may be expressed as 

G t (z*) = (z\U t \0). (24) 

We take the time derivative of the propagator in l|24|l and obtain from the Schrodinger equation l|23[) ihdtGt(z*) — 
(z\H to t(t)Ut\0}- With Hamiltonian JSJl, we arrive at 

d t G t (z*) = -~HG t {z*)+Lz$G t {z*)-itfY,9\e~ iWxt {z\b>V t \Q), (25) 
h A 

with Zf = -iJ2\9\ z \ eluJxt as before. 

As in [24J, we write (z\b\Ut\0) = (z\Utb\(t)\0) with the Heisenberg operator b\(t) — U^ 1 b\Ut- From the corre- 
sponding Heisenberg equation of motion, 

ihdMt) = Ut^Huxifipt = hgle^Lit) (26) 

with the Heisenberg operator L(t) = U i ^ 1 LUt, we find upon integration 

bx(t) = b x -ig* x f dsL( s y^ s . (27) 



We may conclude that 



ft 

(z\b x U t \0) = (z\U t b x (t)\0) = -ig* x I ds e^ s {z\U t L(s)\0), (28) 



where we have used b x 0} = 0. 

Replacing the last term in by expression J2HJ), the evolution equation for the stochastic propagator becomes 



d t G t (z*) = --HG t (z*) + Lz* t G t (z*) -tf J ds a(t - s) (z\U t L(s)\0) , (29) 
with the expected correlation function a(t — s) = J2 X \g\\ 2 e~ lu ' x ( t ~ s \ at zero temperature. 
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Next it turns out appropriate to introduce the operator 

0{t, s) = UtLispf 1 = UtU^LU.Uf 1 (30) 

in the Hilbert space of system and environment. It allows us to express the term under the memory integral in i|29|) 
in the form 

(z\U t L(s)\0) = (z\O(t,s)U t \0). (31) 

Comparing equation (|29[1 with the NMQSD equation (|15|l with replacement <|16|) and recalling equation l|24|) . we see 

that we aim to find an operator 0(t, s, z*) satisfying 

(z\O(t,s)U t \0) = O{t,s,z*){z\U t \0) =0(t,s,z*)G t (z*). (32) 

Relation l|32|) is crucial for finding an evolution equation for 0(t, s, z*) with respect to s, as to be shown in Section llVl 
With expression i|32|) it is clear that equation 129|) is the linear convolutionless non-Markovian stochastic Schrodinger 
equation J5lJ. 

Thus, the operator 0{t,s,z*) in (or JEJ) may be found by investigating the equation of motion for the 
Heisenberg coupling operator 0(t : s) = UtU~ 1 LU s U t ^ 1 in the Hilbert space of system and environment. Depending 
on the application, this approach may offer a more transparent way to establish an expression for the 0(t, s,z*)- 
operator of a convolutionless formulation than the evolution equation 119fl . In particular, it allows us to derive 
evolution equations for 0(t, s, z*) with respect to s, rather than t. 



III. MASTER EQUATION 



Our stochastic approach to open quantum systems will be employed to derive the master equation for the ensemble 
evolution [2^, [H(j . A particularly simple and straightforward route to the master equation is possible as soon as the 
operator 0(t, s, z*) in 116(1 (or 1(32(1 ~) and thus the integrated operator 0(t, z*) from 1(20(1 turns out to be independent 
of the noise z*. This happens to be true for many interesting cases, as shown below. We stress, however, that for the 
main part of this paper, the Brownian motion model of Section ll VI 0(t, z*) does depend on the noise process z% and 
matters arc more difficult. 



A. Density operator evolution 

According to our construction, the reduced density operator p t is given by the ensemble mean over the solutions of 
the stochastic Schrodinger equation (|21|l . We write pt = A4{Pt} with P t = \ipt(z*)) (ipt(z*)\- Upon taking the time 
derivative in (Q) and employing l(2"T)l . we get 

Pt = - l -[H,p t ]+ LM{z* t P t } - M{d{t lZ *)P t } + M{P t z t }tf -M{P t O\t,z*)}L. (33) 
Apparently, this is still far from being a closed evolution equation for p t . Using a version of Novikov's theorem |5l| , 

M{P t z t } =MHds a(t ~ s)-^P t \ , (34) 

which amounts to a partial integration under the Gaussian probability distribution in (|14|) , we may use the replacement 
of the functional derivative by the operator O in ((T^|l and find M{Ptz t } — M {0(t, z*)P t }. Then equation takes 
the form 

Pt = -~[H, p t ] + [L, M{P t O\t, z*)}] + [M{0(t, z*)P t }, L\ (35) 

n 

A convolutionless evolution equation for the non-Markovian dynamics of the reduced density operator may be derived 
from the knowledge of the operator 0(t,s,z*) in 1)16(1 . For instance, as soon as the exact 0(t, s,z*) is independent 
of the noise z\ or as soon as a noise-dependent 0(t, s, z*) can be approximated well by a noise-independent operator 
0(t,s), we find 0(t,s,z*) ~ 0(t,s) and therefore 0(t,z*) — 0(t). In these cases, equation (|35|l is a convolutionless 
closed master equation for p t : due to M{0(t)P t } = 0(t)M{Pt} = 0(t)p t , we find 

Pt - ~[H,pt] + [L,ptOHt)} + [d(t)p t ,tf], (36) 
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Let us emphasize that equation i|36|) is not of Lindblad form in general, not even in the long-time limit due to its 
non-Markovian nature. Nevertheless, as long as we deal with exact master equations - as for most of this paper 
- positivity issues do not appear. By construction, despite their apparent non-Lindblad form, these exact non- 
Markovian convolutionlcss master equations represent (completely) positive maps. However, as soon as the master 
equation follows from an approximate replacement of the operator 0(f, s, z*), positivity is a difficult and delicate 
subject well beyond the scope of this paper. For such investigations on positivity (importance of initial slips etc.) the 
interested reader is referred to the literature El El El mil. 

Clearly, the convolutionless master equation (|36|) covers a much wider class of open system dynamics than the 
Lindblad equation ijSJ). The non-Markovian properties are encoded in a finite width correlation function a(t— s), and so 
in time-dependent coefficients appearing in the master equation. By construction, we already know the corresponding 
stochastic Schrodinger equation l|21|) . Remarkably, here we derive the master equation from the stochastic Schrodinger 
equation. 

We stress that the master equ. (131)1) holds for a noise independent operator 0(f, s,z*) only. Whenever 0(t, s,z*) 
does depend on the noise z* , equ. is still valid yet the step to equ. (fTTfTft fails and matters are considerably more 
difficult. In this case, no general closed evolution equation for pt is known. However, starting from Ij35|l . one may still 
be able to derive an evolution equation for specific cases, as shown in this paper: an important example of this class 
is provided in Section llVI on Brownian motion, where 0(t, s,z*) depends linearly on the stochastic process z* t . 

B. Convolutionless master equations: simple examples 

As noted before, the operator 0(i, s,z*) appearing in (| IL f>|> is generally dependent on the noise z* , so there is no 
general recipe how to derive the corresponding master equation (see our approach in Sec. IV). Before proceeding to 
our main task, in what follows, we review a few examples where either the operator 0(t, s, z*) according to (|16f) may 
be approximated by a noise-independent operator 0(t) or the exact 0(t, s, z*) does not contain the noise z\. In both 
these cases, the corresponding convolutionlcss master equation l|3t)|) is valid. 

1. Weak coupling 

In many interesting cases, the assumption that the interaction between system and environment is weak is a good 
approximation. The action of the functional derivative in may be systematically expanded in powers of the 
interaction H mt H3, |53( • 

oo 

0{t, s ,z*) = Y.9 n On{t,s,z*) (37) 

n=0 

where g represents the coupling strength. To lowest order, one finds 

Mz*) (38) 
and thus in this lowest order the operator 

0(t, s, z*) « O (t, s, z*) = O (t, s) = e -u?(t-«)/*£ e «*(*-«)/» (39) 
independent of the noise. The relevant time-integrated expression reads 

0(f) » (f) = [ ds a(s)e- lHs/h Le tHs ^, (40) 
Jo 

and enters the master equation (|3^)l for the ensemble evolution. Indeed, the resulting equation is nothing but the 
weak coupling so-called "Redfield" master equation 0, |^1| , including an "initial slip" captured by the initial time 
dependence of 0(t). A weak coupling stochastic Schrodinger equation equivalent to (|21|l (with replacement l|4(Jll ') was 
derived in |52( in a formulation that kept the memory integral over the bath correlation function. The convolutionlcss 
formulation (|21|l is easier to handle, as there is no need to store the state vector ip s {z*) at earlier times s < t. 



6Mz*)/Sz: = 



e -iH(t-s)/h Le iH(t- s )/h _|_ _ ^ 
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2. Near Markov 



If the bath correlation function a(t — s) falls off rapidly under the memory integral in l|15(l . an expansion of the 
functional derivative in terms of the time delay (t — s) is sensible [22I I23I Isoj . With A n (t) — ds s n a(s) and 
n = 0, 1, 2, . . ., we find for the relevant integrated operator l|2U|) to first order, 

6(t) n A (t)L + At(t) (-^[H,L] + A (t)[L,tf]Lj , (41) 

neglecting contributions A n (t) with n > 2. Once again, to the order accepted, we see the operator 0(t, s,z*) turns 
out to be independent of the noise z\ and the master equation ll.'UiH immediately applies (see also |50j). 

Of particular interest is the standard Markov limit: the correlation function may be replaced by a delta function, 
a(t — s) — 7<5(i — s), with some constant 7. The only relevant term in 1)41(1 is the Ao-term which may be replaced by 
the constant O — ^L. As noted earlier, the stochastic Schrodinger equation \2l\i (and thus, IjlSI) ) reduces to the linear 
version of the Markov quantum state diffusion equation @ . Correspondingly, the master equation H36|) is nothing but 
Lindblad's equation J3J. The next order Ai-term in expansion l|41[l turns out to be relevant for the high temperature 
limit of the quantum Brownian motion model |23| and quite generally leads to a theory of "post-Markov" evolution 




We remark that very often, Markov and weak coupling approximation are only meaningful in combination, referred 
to as Born-Markov approximation. 



3. Soluble models 



For some choices of system Hamiltonian H and coupling operator L in the interaction Hamiltonian © , the functional 
derivative may be replaced by an operator 0{t, s, z*) without any approximation |18|. Two examples are mentioned, 
where the O-operator turns out independent of the noise z* . The third, more complicated case of Brownian motion 
of a harmonic oscillator will be the main subject of the following Sections. 

The first simple example is an harmonic oscillator with H — fla^a, coupled to the environment through a rotating- 
wave-type coupling, L = a in the total Hamiltonian J5J|. It turns out |18| that the functional derivative in l|16|l may 
be replaced by an operator 0(t, s) being proportional to the annihilation operator without any approximation, 

0(t, S ) = ^-a. (42) 
c(t) 

Here, c(t) is a complex function satisfying the equation of motion 

c(t) + iVLc{t) + [ ds a(t - s)c(s) = 0, (43) 
Jo 

involving the bath correlation function. We recognize the damped motion of the amplitude of the oscillator. The 
integrated operator 0(t) thus becomes 

0(t) = C(t)a (44) 

without any approximation, where C(t) — f Q ds c(s)a(t — s) / c(t) . The resulting exact convolutionless master equation 
[f34| thus follows from the general expression (|36ll : 

pt = -i[na^a,pt}+C*(t)[a,pt^]+C(t)[ap t ,a^] (45) 
= -*[(n + lm{C(t)}) a^a, p t ] + Re{C(t)} ([a, p t a^} + [ap t ,o)]) 



A similar result may be derived for a damped two- level system with Hamiltonian H = Qo~ z /2 and the rotating- wave 
type coupling L = er_ [l8| . Similar to the case of the damped (rotating- wave) harmonic oscillator above, one finds 

0(t, S ) = ^L (46) 
c(t) 



with the very same function c(t) from l|43|l . Thus, an exact convolutionless master equation (|36|) of a form similar to 
(|45|l follows, with a replaced by cr_ in the non-unitary part of the evolution equation. 
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To summerize, in all the cases presented in this subsection it is possible to replace the functional derivative in 1|16|) 
by an operator 0(t, s) that is independent of the noise z t *. In such cases, the ensemble mean follows the convolutionless 
master equation (I36|) . In the weak coupling case with 0{t) from (|39(l . for instance, we recover the so-called Redfield 
master equation [47], 03 including an initial slip. In the Markov case, keeping only the lowest order in expansion l|4ip. 
we recover Lindblad's equation Q for the ensemble evolution. For the two exactly soluble models with 0(t, s) from 
(|42JI or H46|l . equation (|36|l is the exact convolutionless master equation. In these latter cases, asymptotically for large 
times t, when the time dependent function C(t) approaches a constant, the exact master equation turns into one of 
the Lindblad class (0- In the following Sections we discuss yet another exactly soluble case, where the convolutionless 
master equation is more involved due to the dependence of the operator 0(t, s, z*) in l|16|) on the noise z t *. 



IV. QUANTUM BROWNIAN MOTION 
A. Brownian motion model 



Brownian motion of a harmonically bound particle may be obtained from the Hamiltonian 



A 



of an oscillator with mass M and f req uen cy i l , coupl ed to an environment of harmonic oscillators through its position 
q [22l 0, EE EE 03 EE HE EE EE EE EE EH EE 03 • In this form, the usual "counterterm" arising from the coupling 
should be understood to be included in the harmonic potential, see |E for more details. 

The Brownian motion Hamiltonian (|47fl is of the form JSJ and thus may be treated with our stochastic approach 

to open quantum systems. With q\ = (b\ + b\j we identify the quantities entering the basic Hamiltonian (JSJ) 



to be g\ = G\\l 2m ^ and the coupling operator is L = = q/h. The system Hamiltonian is the harmonic 



v 2 ! 

H= — + -Mn 2 q 2 . (48) 
2M 2 K ' 

For the quantum bath correlation function JJjJ of this model at temperature T one finds the force-force correlation 



a(t - s) = ((B(t) + B\t))(B(s) + Bt( s ))) env (49) 

2k B T, 



~^ 2m\LO\ 



coth I — — ^— I cos u)\(t — s) — ism.ui\(t — s) 



Introducing a spectral density of bath oscillators 

V 

2m\uj\ 



A 

one usually writes 



/>oo 

a(t — s) = h I dw J 
Jo 



coth | | cosluU — s) — i smujtt — s) 



(51) 



Often the so-called Ohmic case is chosen with J(ui) = M^ujf c {u>/K) with some cutoff function f c (x), e.g. f c (x) ~ e 
and A a cutoff frequency. We stress, however, that the following results are valid for any spectral density J(u>). 



B. Master equation 

It is well known that the model l|47|) allows the derivation of an exact convolutionless master equation for the 
reduced density operator |34l l36l l37| . It may be written in the form 



Pt = ±[H,p t ] + ^ll 2 ,Pt] + |f[9, {p,Pt}} + b*]] - k,Pt}] IWI 
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with real time dependent coefficients a(t) , b(i) , c(t) and d(t). The physical meaning of these terms (drift terms a(t) 
(frequency renormalization), bit) (damping term) and diffusion terms c(t),d(t)) become s appar ent from the Wigner 
representation of equation (j52(l , for which the reader is referred to the literature 0, 1371 l39t I4II Ep . 

In the derivation of l|52[l was based on the generator of the time evolution, while in |37Ll3q. the authors used 
path integrals. In Section we show how the master equation ll52")) follows directly from the stochastic Schrodinger 
equation which we discuss next. 



C. Stochastic Schrodinger equation 

Our stochastic Schrodinger equation approach allows a rigorous treatment of model 14711 . The replacement of the 
functional derivative with an operator 0(t, s, z*) in (|16|) can be established without any approximation, and thus a 
convolutionless exact stochastic Schrodinger equation (|21|l may be found. Here, however, we have to allow for an 
explicit noise dependence. It turns out [l]| that the ansatz 

0(t, s, z*) = i (f(t, s)q + ^g(t, s)p - ^ J ds'j(t, s, s')z* s ,) , (53) 

with complex functions fit, s), g(t, s), and jit, s, s') to be determined is a solution of the general evolution equation 
(|19|l . We learn from this example that in general, the functional derivative in l|16|l may indeed introduce a dependence 
on the whole history of the noise z% . 

The functions fit, s), g(t, s) and jit, s, s') in l|53f) have to satisfy the evolution equations 

d t f(t,s) = n g (t,s)-2ig(t,s)F(t)+if(t,s)G(t)+iJ{t,s), (54) 
d t g(t,s) = -Qf(t,s)-ig(t,s)G(t), and 
d t j(t,s,s') = -ig(t,s)J(t,s'), 

where we introduced the integrated functions 

F{t) = ikh f dsa ^- s )f^ s )' (55) 
1 /■< 

G(t) = j-j^ / dsa{t-s)g{t,s), and 
1 /■* 

J(t,s) = jjQfo J dsa(t-s)j(t,s,s). 

While /, g, and j are dimensionless complex functions, the integrated expressions F, G, and J are defined such as to 
have dimension inverse time: typically, these latter functions turn out to be proportional to a damping rate 7. The 
evolution equations (|54|l have to be solved with boundary conditions 

f(t = s, s) = 1, g(t = s, s) = 0, j(t = s, s, s') = 0, and (56) 
j(t = s,s',s) = -g(s,s'), 

for all s and s' . The first three conditions arise from the initial condition (|18|l for the operator 0(t = s, s, z*) = L = 
q/Ti, the last condition arising from the solution of l|19|) . The relevant integrated operator (|20|) turns out to be 

0(t, z*) = MQF(t)q + G(t)p - i [ dsJ{t, s)z* s (57) 

Jo 

with time dependent coefficients from (|55H . 

The resulting linear, convolutionless non-Markovian stochastic Schrodinger equation (|21() for the Brownian motion 
of a harmonic oscillator (|48|l with a coupling to the environment through position L = q/h thus reads 

hdtipt = -iHip t + q(zt -0(t,z*))i> t (58) 
= -iHipt + qiz^ -MClF(t)q-G(t)p + i [ ds J(t, s)z*] ip t . 



The process z t * is a complex Gaussian process with the finite temperature correlation function a{t — s) from (|51|l . Just 
like the exact master equation l|52|) . the exact stochastic NMQSD equation (|58|l is valid for arbitrary bath correlation 
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function a(t — s), and thus for any temperature, environmental spectral density, or coupling strength between the 
harmonic oscillator and its environment. 

This example demonstrates clearly that the evolution of quantum trajectories ipt(z*) may indeed depend on the 
whole history of the noise. Simulations of the Brownian motion model discussed here may be found in [55j. 

From our construction, it is clear that the ensemble of quantum trajectories of the exact equation l|58l) evolves 
according to the exact master equation (59 ■ I n the next Section we set out to establish this connection directly. 

As 0(t, z*) depends on the noise, the simple result (|36[) for the evolution of the reduced density operator p t is not 
applicable, and further effort is required. As a first step, we show in the next Section how a Heisenberg operator 
approach helps to gain more insight into result l|53[l for the operator 0(t, s, z*) replacing the functional derivative. 



D. Heisenberg method 

Major clarification is achieved once we investigate the dependence of the functions f(t,s),g(t,s), and j(t,s,s') in 
(|53|l on s rather than t. The key step arises from the Heisenberg operator approach as explained in Section III Dl see 
also |2^. Here, the challenge is to find a suitable expression for 

Q(s) = (z\U t q(s)\0) = hO(t,s,z*){z\U t \0), (59) 

similar to expressions (|31H and i|32|l . Clearly, Q(s) depends on the time t, but we here regard it as a function of s and 
thus may investigate the s-dependence of the desired operator 0(t, s, z*). Not surprisingly, the Heisenberg equation 
of motion for q(s) leads to d s Q(s) — (z\Utd s q(s)\0) = (z\U t p(s)\0) and therefore we investigate the second order 
derivative. Note that required final values are 



Q(s = t) = q{z\U t \0) (60) 
1 

M' 



dsQ(s)\ s=t = YjPizlUtlO). 



For the second order derivative we find from (|47(l with g\ = w 2m ^ A 

3 s 2 Q(s) = -n 2 Q(s) -j-^gxe-^ s (z\U t b x (s)\0) - £ g x e^ s (z\U t b{(s)\0). (61) 

A A 

As before in Section (f II D(l . integrating the Heisenberg equation of motion for the environmental annihilation 
operator b\ (s), we see that 

(z\U t a x (s)\Q) = -^g x J° ds' e^ s '(z\U t q(s')\0), (62) 

where use is made of the fact that the initial b\(0)\0) — b\\Q) = 0. 

The adjoint b x (s) in (|61|) is dealt with in similar fashion. Here, however, we make use of the fact that at the final 
time t the overall expression becomes simple: (z\U t b x (t)\0) = z* x {z\U t \Q)- Thus, we integrate the Heisenberg equation 
of motion for b\{s) given the final value at s = t to get b\(s) = b x (t) — j-g\ J ds' e~ tuJxS q{s'), and find 

(z\U t b\(s)\0) = z* x (z\U t \Q) - l -g x J ds' e^' (z\U t q(s')\0) . (63) 

All that remains to be done is combining the results l|61(l , (162t , and l|63|l to obtain a linear second order differential 
equation for Q(s): 

8 2 s Q{s) + n 2 Q{s) (64) 

~WhJ ds ' a<yS ~ S '^ S ') ~ ~WhJ s ds ' a<yS ' ~ s ^ s ') = -jj z s( z \ Ut \°^ 

for ,s G [0, t) with a fixed final time t, and final values (|60|l . 

With expression l|59|l . the derived second order differential equation for Q(s) translates into an equation for the 
desired operator 0(t, s, z*). Considered as a function of the time s, at fixed t, we find 

d 2 s O(t,s,z*) + n 2 0(t,s,z*) (65) 

-m f ds ' a{s - s ' )0 ^ s ' ^ - m f s ds ' a{s ' - s)0(< ' - ~ik z ° 
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with final values from l|60|) 



0(t,a=t,z*) = ±q, and d s O(t, s, z*)\ s=t = ^p. ((Hi! 



As before in (|53[l . we write the operator 0{t, s, z*) in the form 

0(t,s,z*) = ^ (<t>t(s)q + J^M S )P~ J o ds'xt{s,s')z* s ,^j , (67) 

where we introduced three s-dependent complex functions 4>t(s), ipt{s) and Xt( s ; s ') at a fixed time i. Obviously, from 
I|fi5jl. the first two of them satisfy the homogeneous equation 

d 2 s 4> t { S ) + tfus) -jfi-f ds " «( s - *")M*") ~wif s ds " a(s " - s) ^ (s,,) = °- (68) 

While for 4> t we demand the final values 

<p t (s = t) = l and 0»&(s)| a=t = O (69) 

to satisfy Ij66(l . ipt{s) has final values 

Ms = t) = and S! B ^(*)L =t =fi. (70) 
The function Xt(s, s/ ) in (|67|l is nothing but Green's function that satisfies 

d 2 sX t(s,s') + n 2 Xt ( s , s ') (7i) 

' f S ds" a{s - s") Xt (s", s') ~ TTT / ds " a ( s " ~ s )Xt{s" , a') =SlS(s- s') 



Mh J Q v wv ' ' Mh^ 

with final values 

Xt (s = t,s')=0 and d sX t(s,s% =t = (72) 

for all s' e [0,t]. 

The result (|67|l for the operator 0(t, s, z*) reflects our ansatz ()53|l that solved evolution equation (|19|l with respect 
to t. Clearly, we identify the relations f(t 7 s) = (f>t{s), g{t,s) — ipt{s), and j(t,s,s') = xt(s,s'). All conditions f5fi|l 
on the functions f,g,j are met by (|69|) . I|7U|I . and (|72|l . One may even derive their evolution equations l|54|l regarded 
as a function of t from the corresponding evolution equations (I68|l and (|71|l . when considered as functions of s, see 
appendix^ for an example. For numerical applications as in [55) . the evolution equations with respect to t are the 
ones of interest. In order to derive the convolutionless master equation for this quantum Brownian motion model, 
however, it turns out that the evolution equation of 0(i, s, z*) with respect to s from l|65|l is the right starting point, 
as will be explained in the next section. 

V. QBM MASTER EQUATION 

In the last Section we have established the evolution equation with respect to s for the operator 0(t, s, z*) entering 
the exact convolutionless stochastic Schrodinger equation l|58|) . We are now in the position to derive the corresponding 
convolutionless master equation from the general expression (|35ll . 

dtPt = ~[H,pt] + \[q,M{P t O\t,z*)}\ + ±[M {0(t,z*)P t },q\ (73) 

which for the Quantum Brownian motion case we display here with coupling operator L = q/Ti. 

We saw previously how knowledge about the operator 0(t, s, z*) replacing the functional derivative in our stochastic 
Schrodinger equation, puts us in a position to derive a closed convolutionless master equation. In particular, as soon as 
the operator 0(t, s, z*) is independent of the noise, we get A4{0(t)P t } — 0(t)p t and (|73|l provides the convolutionless 
master equation (|36() for p t . In our current example, however, 0(t, z*) does depend on the noise and we show next 
how to proceed in this case. 
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Starting point is the evolution equation i|65|l for the operator 0(t, s, z*) as a function of s. For the master equation 
(|73[) . we need an expression for the ensemble mean 

M {6(t, z*)P t } = [ dsa(t- s)M{0{t,s,z*)P t } = [ dsa(t - s)R(t,s) (74) 
Jo Jo 

with 

R(t,s) = M{0{t,s,z*)P t }. (75) 

Upon taking the mean Ai{. . . Pt] of equation I|tj5fl . we see from its definition (|75|) that R(t,s) satisfies the second 
order differential equation 

d 2 R{t, s) + n 2 R(t, s) (76) 

~m[ ds ' < s ~ w> *'> ~ m [ ds ' a{s ' ~ s)Rit > s>) = ~m M ^ } - 

Next we have to find an expression for A4{z*P t }, which again follows from Novikov's theorem [5l| as in 1)34(1. We 
may express the resulting functional derivative in terms of the operator 0(t, s, z*) and get 

M{z*P t }= f ds'a*(s-s')M{P t O*(t,s',z*)} = f ds'a*(s - s')R*(t,s'). (77) 
Jo Jo 

A closed equation for the operator R(t, s) as a function of s results: 

d 2 R(t,s) + n 2 R(t,s)--^- [ ds' a{s-s')R(t,s') (78) 



Required final values are 



Mh J 



^ ds' a(s' -s)R(t,s') = -^r I ds>a*(s-s')RHt,s>). 



R(t,s = t) = -qp t , and d s R{t,s)\ s=t = -jj^PPt, (79) 



as one can easily see from (16611 . 

In fact, it is clear that equation (|78|l with final values (|79|l may be satisfied by an expression of the form 

R(t, s) = ^ jfc(i, s)qp t + s )PPt + m (t, s)p t q + J^n(t, s)p t p^ , (80) 

where k(t, s), £(t, s), m(t, s), and n(t,s) are complex functions whose s-dependence follows from equation (|78|l with 
appropriate final values at s = t from (|79|l . As an example, we get k(t, s = t) = 1 and d s k(t, s)\ s=t = 0. In Appendix 
[5]we show how to determine all four functions. 

Now that we know R(t, s) = Ai{0(t, s, z*)P t } from (|80|l . the closed master equation for the reduced density operator 
follows immediately from the general expression (|73|l with (|74l) . After some rearrangements, we indeed arrive at the 
more familiar form l|52|) . with time dependent coefficients 



2 



a(t) = -Im ( / dsa(t- s)(k(t, s) + m(i, s)) ) (81) 



2 < '* 



c(t) = ]^ Rc (^ dsa(t-s)(n(t,s)-i(t,s)) 
d(t) = Re (^J ds a(t - s)(k(t, s) - m(t, s)) 



In appendix [B] we discuss these expressions in more detail, establishing the connection to earlier derivations of the 
master equation l|52|) . Crucially, we here show explicitely that it follows from the corresponding convolutionless 
stochastic Schrodinger equation (J5SJ. 
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VI. CONCLUSION 



Exactly soluble models are valuable tools. They allow us to discuss approximation schemes, show the transitions 
in the qualitative behaviour when parameters are changed over a wide range, and last not least provide insight into 
different theoretical approaches. 

Our non-Markovian stochastic Schrodinger equation approach offers a new method for handling open quantum 
systems. In this paper we use this framework as a theoretical tool to derive convolutionless non-Markovian master 
equations. Most interestingly, the exact master equation for quantum Brownian motion of a harmonic oscillator, 
coupled to a bath of oscillators is derived from the corresponding convolutionless non-Markovian stochastic Schrodinger 
equation. In both approaches, non-Markovian properties are encoded in time-dependent coefficients. 

The problem how to describe non-Markovian open quantum system dynamics efficiently is a difficult one. As further 
underlined in this papers, wc believe that our non-Markovian stochastic Schrodinger equation approach offers a useful 
alternative framework to the existing ones. 
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APPENDIX A: EVOLUTION EQUATION WITH RESPECT TO TIME t 



It is interesting to see the connection between the evolution equations l|54|) with respect to t for f(t, s) = 
4>t(s), g(t, s) — -0t( s ) an d 5, s') = Xt{s,s') and the differential equations lj6l3|) and l(7T)l for the same functions 
with respect to s. In fact, here we show how to derive the former from the latter. 

Consider g(t,s) = ipt(s) as an example. We showed in Section TlV Dl that this function solves I|ti8[) which we here 
display again for convenience: 

d 2 s Ms) + u?Ms) -jrh[ ds " a{s ~ s ")^( s ") -ml ds " a[s " ~ s)Ms " ] = °' (A1) 

Taking the time derivative of this equation with respect to t, and using the final condition ip t (s = t) = from H70JI . 
we find that dtipt(s) satisfies the very same equation (|A1|) . Therefore, as {0t(s), V't(s)} forms a basis of solutions, 
there exists a relation of the form 



d t Ms) = ci (t)&(«) + c 2 (t)Ms) 



(A2) 



with suitable coefficients ci(t) and C2(t). These have to be determined from the final values at s = t. From l|t)9|) and 
ffTUfl we get 



ci(i) = d t i>t(s = t) 



(A3) 



These expressions are easily evaluated: first, from ipt(s = t) = 0, we get dtipt(s = t) + d s ipt(s)\ s=t = and thus 
c i{t) = — d s tjjt(s)\ s=t = — O according to {7UJ. Secondly, we use the trivial identity 



d„Ms) = Ws)| s=4 - / ds" (9 s 2 Vt(s)L s „) 



(A4) 



n 



ds" (d 2 s Ms)\ 
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Take the derivative with respect to t and then set s = t to find 

= -d 2 s Ms)\ s=t (A5) 

= n 2 M* = t )-j^Jl ds " «( f ^ 

= -j^rj* ds" a(t-s")Ms'% 

where use is made of the evolution equation (|A1I) and the condition l|70() . Thus, with (|A3fl . C2(t) = — M ^ ft J ds" a(t — 
s")ip t (s") which may be written as C2(t) = —iG{t) using definition l|55|l . We see that with C\(t) — — Cl and C2(t) — 
—iG(t), relation 1A2(I is nothing but the evolution equation for ipt{s) = g{t, s) in i|54J) as derived from the "consistency 
condition" (|T7jl . 



APPENDIX B: TIME DEPENDENT COEFFICIENTS 



From the second order differential equation (j78(l for the operator R(t, s) which we here display again for convenience, 
R(t, s) = i jfc(£, s)qp t + s )PPt + s)p t q + s )PtP j , (Bl) 

we get (coupled) differential equations for the coefficients k(t, s), £(t, s), m(t, s), and n(i, s) as functions of s. The final 
values at s = t for all four functions are determined from the conditions (|79f) on R{t 1 s). 

It turns out useful to replace the four unknown functions by the linear combinations appearing in the final result 
(|8ip. which we define as 

w(t,s) = k(t,s) +m(t,s) (B2) 
x(t,s) = l(t, s) + n(t, s) 
y(t,s) = k(t, s) — m(t, s) 
z(t,s) = n(t, s) — £(t, s) 

They all satisfy uncoupled second order differential equations. For w(t, s) we find from Ij78(l 

i f s 

d 2 s w(t,s) +n 2 w(t,s) / ds'a(s- s')w(t,s') 

Mh J 

' ' a{s' -s)w(t,s') = -/ ds'a(s' - s)w*{s'), (B3) 



MhJ s v ' v ' Mh_ 

and the identical equation for x{t, s). Required final values are 

w(t,s = t) = l, d s w(t,s)\ s=t = (B4) 
x(t,s = t) = 0, d s x(t, s)| s=t = Cl. 

For y(t, s) (and also z(t, s)) we get a similar equation with a different sign on the right hand side: 

d 2 s y(t, s) + Vl 2 y(t, s)-^£ ds'a(s - s')y(t, s') 

~ikf s a{s ' ~ s)y{t ' s>) = mSl ds ' a{s ' ~ s)v * {s% (B5) 

here with final values 

y(t,s = t) = l, d s y(t,s)\ s=t =0 (B6) 
z(t,s = t) = 0, d s z(t,s)\ s=t = -Q. 

It turns out that matters simplify once we separate real- and imaginary parts of the four functions w,x,y, z. For 
the bath correlation function we write 

a(t-s) =v(t-s)+ihr](t-s). (B7) 
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Thus, microscopically, from its definition @ we have 

r)(t — s) = — } — smu)\(t — s) = — / duJ(uj) sinw(i — s). (B8) 

2m\cj\ J 

This kernel turns out to be crucial for the classical equation of motion of the underlying model l|47|l . as will be clarified 
shortly. 

For the imaginary part of w(t, s) we get from (|B3|I 



2 /"* 

d 2 s w I (t,s) + n 2 w I (t,s) + — ds'r)(s-s')w 1 (t,s')=0 (B9) 

M J s 

and the very same equation for the imaginary part of x(t,s). With final values wi(t,t) = 0, d s Wi(t, s)| s =j = and 
xj(t, t) = 0, d s xj(t 7 s)| s= t = the solutions are trivial: 

Wi(t,s) = for all s (BIO) 
Xi(t, s) = for all s. 



Next we consider the real parts Wait, s),x R {t, s). With (|B10(I and (|B3J) we arrive at 



d 2 s w R {t,s) + n 2 w R (t,s) + ^- f ds'ri(s-s')w R (t,s')=0 (Bll) 

and the same equation for x R (t, s). Required final values are 

w R {t,s = t) = l, d s w R (t,s)\ s=t = (B12) 
x R (t,s = t) = 0, d s x R {t,s)\ s=t = n. 

Similarly, we find for the real and imaginary part of the function y(t, s) the following equations: 

d 2 y R (t,s) + n 2 y R (t,s) + |- / d8'ri(8'-8)y R (t,s') = (B13) 



2 r „ 2 '■' 



M 

d 2 s y I (t,s) + n 2 y I (t,s) + j^ J ds'r)(s-s') yi (t,s') = J ds'v(s - s')y R (t, s'). 
As before, real and imaginary part of z(t, s) satisfy the very same equations. Required "final values" are 

(B14) 



VR{t,t) 


= 1, 


dsyn(t,s)\ s =t = 


yi{t,t) 


= 0, 


d s yi(t,s)\ s=t = 


ZR(t,t) 


-o, 


d s z R (t, s)\ s=t = -SI 


Zl{t,t) 


= 0, 


d s zi(t, s)\ s=t = 0. 



The real parts of all four functions may nicely be expressed using the special solution q(s) of equation I|B11|) 

d 2 q{s) + n 2 q(s) + A / ds'r)(s - s')q(s') = (B15) 



with initial values 

g(0) = 0, and d s q(s)\ s=0 =a (B16) 

Equation (|B15|I is nothing but the classical equation of motion for the position q(s) of the underlying model Ij47|l . 

provided the environmental oscillators all start with initial values q\(0) = 0, q\(0) = 0. Note that q(s) satisfies the 
very same equation l|B15|) with q(0) = 57 and q(0) = 0. 

It is customary to write equation (|B15(I in the more familiar form 

d 2 q(s) + n 2 q{s) + [ ds' 7 cl(s - s')q(s') = (B17) 
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with the classical damping kernel 7ci(s — s') defined through jjrj(t — s) = dt"f c \(t — s). 
In terms of this classical solution q(s) simple inspection shows that 



- miftmr - (B18) 

*■•>-■«<«■•> - a W^- 

yn(t,s) = —q(t-s), 
zn(t,s) = q(t - s). 

While the imaginary parts of w(t, s) and x(t, s) are zero, the expressions for yi(t, s) and zj(t, s) are more involved. 
We find 

Vi(t, s) = f ds> f ds"u(s' ~ s")q(t s")q(s - s') (B19) 



2 



ds' / ds"v(s' -s")q{s")[w(t,s)q(s')+x{t,s)q(s')/fi\, 



MVL 2 fi 

Zl(t,a) = [ r/,' / -.")</(/ -s"), l{ , .s') 

2 ' ds' I ds"u{s' -s")q(s")[w(t,s)q(s') + x(t,s)q(s')/n]. 



urn 

The first line in each expression is simply the inhomogeneity on the right hand side of equations l|B13(l , in convolution 
with Green's function q(s — s')/tt. The second lines, respectively, are solutions of the homogeneous equations and 
serve to satisfy the final conditions (|B14(I . 

Finally, let us make contact to earlier derivations of the time dependent coefficients. For brevity we concentrate on 
the drift coefficients a(t) and b(t) in the master equation As w(t,s) and x(t, s) are real, using l|gT|) and (|B7|I . 

a(t) and b(t) may be written in the form 

a(t) = 2 / dsrj(t - s)w R (t,s) (B20) 



With expressions 

jnilt , and using the evolution equation (|B 1 5|l for both, q(s) and q(s), we may write the drift 
coefficients in the form 

a(t) = -A^ 2 + M«^i»l, (B21) 



b(t) 



q(t)q(t) - q{t)q\t) 



Expressions similar to these may also be found in the literature [33. l37l 13^. lioj . 



[1] C. H. van der Wal et. al. Science, 290, 773 (2000); J. R. Friedman et. al. Nature 406, 43 (2000); 

[2] M. Arndt et. al Nature 401, 680-682 (1999); L. Hackermiiller et. al. Phys. Rev. Lett. 91, 90408 (2003). 

[3] F. Schmidt-Kaler et. al. Nature 422, 408 (2003); D. Leibfried et. al., Nature 422, 412 (2003). 

[4] H. Carmichael, An open system approach to Quantum Optics, (Springer, Berlin) (1993). 

[5] C. W. Gardiner and P. Zoller Quantum Noise, 2nd ed. (Springer, Berlin) (2000). 

[6] U. Weiss, Quantum dissipative systems, (World Scientific, Singapore), 2nd edition, (2000). 

[7] H. Carmichael, Statistical methods in quantum optics I, (Springer, Berlin) (1999). 

[8] H.-P. Breuer and F. Petruccione The theory of open quantum systems, (Oxford University Press, Oxford) (2002). 
[9] J. Dalibard, Y. Castin, and K. M0lmer, Phys. Rev. Lett. 68 580 (1992). 



18 



N. Gisin and I. Percival, J. Phys. A 25, 5677 (1992); 26, 2233 (1993). 
C. Gardiner, A. Parkins, and P. Zoller, Phys. Rev. A 46 4363 (1992). 
H. M. Wiseman and G. J. Milburn, Phys. Rev. A 47, 1652 (1993). 
M. Plenio and P. L. Knight, Rev. Mod. Phys. 70, 101 (1998). 

G. Lindblad, Comm. Math. Phys. 48, 119 (1976). 

V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, J. Math. Phys. 17, 821 (1976). 
W. T. Strunz, Phys Lett. A 224, 25 (1996). 
L. Diosi and W. T. Strunz, Phys. Lett. A 235, 569 (1997). 
L. Diosi, N. Gisin and W. T. Strunz, Phys. Rev. A 58, 1699 (1998). 
W.T. Strunz, L. Diosi, and N. Gisin, Phys. Rev. Lett. 82, 1801 (1999) 
M. Jack, M. Collet, and D. Walls. Phys. Rev. A 59, 2308 (1999). 

H. Breuer, B. Kappler, and F. Petruccione, Phys. Rev. A 59 1633 (1999). 
T. Yu, L. Diosi, N. Gisin and W. T. Strunz, Phys. Rev. A 60, 91 (1999). 
W. Strunz, L. Diosi, N. Gisin and T. Yu, Phys. Rev. Lett. 83, 4909 (1999). 
J. Cresser, Laser Phys. 10, 1 (2000). 

J. Gambetta and H. M. Wiseman, Phys. Rev. A 66, 012108 (2002). 
S. Nakajima, Prog. Theor. Phys. 20, 948 (1958). 
R. Zwanzig, J. Chem. Phys. 33, 1338 (1960). 

F. Schwabl and W. Thirring, Ergeb. exakt. Naturwiss. 36, 219 (1964). 

N. G. van Kampen, Stochastic Processes in Physics and Chemistry, 2nd ed. (North-Holland, Amsterdam), (1992). 
H.-P. Breuer, B. Kappeler, and F. Petruccione, Phys. Rev. A 59, 1633-1643 (1999). 
R. Feynman and F. L. Vernon, Ann. Phys. (N.Y.) 24, 118 (1963). 
P. Ullersma, Physica 32, 27 (1966). 

A. O. Caldeira and A. J. Leggett, Physica A 121, 587 (1983); 

F. Haake and R. Reibold, Phys. Rev. A 32, 2462 (1985). 
H. Grabert, P. Schramm, and G.-L. Ingold, Phys. Rep. 168, 115 (1988). 
W. G. Unruh and W. Zurek, Phys.q Rev. D 40, 1071 (1989). 

B. L. Hu, J. P. Paz and Y. Zhang, Phys. Rev. D 45, 2843 (1992). 
B. L. Hu, J. P. Paz and Y. Zhang, Phys. Rev. D 47, 1567 (1993). 
J.J. Halliwell and T. Yu, Phys. Rev. D 53, 2012 (1996). 

G. W. Ford and R. F. O'Connell, Phys. Rev. D 64, 105020 (2001). 
L. Diosi, Physica A 199, 517 (1993); Europhys. Lett. 22, 1 (1993). 
E. Calzetta and A. Roura, and E. Verdaguer, Physica A 319, 188 (2003). 
W.T. Strunz and F. Haake, Phys. Rev. A 67, 022102 (2003). 

S. John and T. Quang, Phys. Rev. Lett. 74, 3419 (1995); B. Garraway, Phys. Rev. A 55, 4636 (1997); N. Vats and S. 
John, Phys. Rev. A 58, 4168-4185 (1998); M. W. Jack, M. J. Collett, and D. F. Walls, Phys. Rev. A 59, 2306 (1999); G. 
M. Nikolopoulos and P. Lambropoulos, Phys. Rev. A 61 053812 (2000); M. W. Jack and J. J. Hope, Phys. Rev. A 63, 
043803 (2001). 

[45] J. J. Hope, Phys. Rev. A 55, R2531-R2534 (1997); H.-P. Breuer, D. Faller, B. Kappeler, and F. Petruccione, Phys. Rev A 
60, 3188-3196 (1999); J. J. Hope, G. M. Moy, M. J. Collett, and C. M. Savage, Phys. Rev. A 61, 023603 (2000). 
S. Gnutzmann and F. Haake. Z. Phys. B 101, 263 (1996). 
P. Gaspard and M. Nagaoka, J. Chem. Phys. Ill, 5676 (1999). 
V. Bargmann, Commun. Pure Appl. Math. 14, 187 (1961). 
T. Yu (to be published). 

T. Yu, L. Diosi, N. Gisin and W. T. Strunz, Phys. Lett. A 265, 331 (2000). 
A. Novikov, Sov. Phys. JETP 20, 1290 (1965). 
P. Gaspard and M. Nagaoka, J. Chem. Phys. Ill, 5668 (1999). 
W. Pfersich and W. T. Strunz (to be published). 
R. R. Puri and S. V. Lawande, Phys. Lett. A 62, 143 (1977). 
W. T. Strunz, Chem. Phys. 268 237 (2001). 



